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We have conducted temperature and frequency dependent transport measurements in amorphous 
NbjjSii-a; samples in the insulating regime. We find a temperature dependent dc conductivity 
consistent with variable range hopping in a Coulomb glass. The frequency dependent response in 
the millimeter-wave frequency range can be described by the expression cr{ijj) oc {—luj)" with the 
exponent somewhat smaller than one. Our ac results are not consistent with extant theories for the 
hopping transport. 
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The study of conducting processes in disordered insu- 
lators has played an important role in the history of the 
physics of disordered materials yet retains considerable 
intrinsic interest. These materials are sometimes called 
"Coulomb glasses" for the presumed complexity of the 
interactions between localized electrons. Efforts to un- 
derstand the conducting mechanisms in the insulating 
phase continue as a variety of mechanisms are possible — 
including thermal activation to a mobility edge, thermal 
activation to a neighboring localized state, and variable 
range hopping. For materials well within the insulat- 
ing phase, variable range hopping is likely to dominate 
the transport; when a delocalization transition is ap- 
proached, however, it is not obvious which process will 
dominate In addition, the role of strong electron- 
electron interactions, which is known to be important in 
the vicinity of the metal-insulator transition has 
not been fully explored in the insulating state. 

Theories based on the variable range hopping mech- 
anism either without strong effects due to electron- 
electron interactions (Mott hopping f^) or with them 
(Efros-Shklovskii (E-S) hopping @,§) are well estab- 
lished in explaining various aspects of the dc and low 
frequency dielectric response of lightly doped semicon- 
ductors and very disordered alloys. They each produce 
characteristic temperature dependences: 

Mott : adc oc exp[-(To/r)i/('^+i)] (1) 
E-S : arfe(xexp[-(ro/T)i/2]. (2) 

Here d is the effective spatial dimension for the hopping 
electrons and Tq is a characteristic temperature scale that 
is related to the level spacing in a volume (where ^ is 
the localization length) in the Mott theory and to the 
Coulomb interaction energy scale in the Efros-Shklovskii 
theory: fc^To = b^/k^ (where k is the dielectric con- 
stant). Note that these expressions are valid for T <C Tq. 



Frequency dependent conductivity experiments are 
particularly useful in distinguishing between different 
conducting mechanisms since they directly probe the 
quantum state of the electrons (especially in the high 
frequency limit where Hlo > kgT). The frequency- 
dependent conductivity has been measured in the disor- 
dered insulating state from audio to microwave frequen- 
cies [^,|o|, far- infrared frequencies ||ll|,|l^, and beyond 
JlSf . For investigating the quantum behavior of the hop- 
ping conductivity, however, measurements at frequencies 
lower than far-infrared (below the frequency kBTo/h or 
the range where other excitations such as phonons are 
important) but still high enough to be in the quantum 
limit huj > k^T, i.e., low temperature microwave and 
millimeter-wave conductivity experiments, are most rel- 
evant. Here we report results for the complex conduc- 
tivity in this frequency range at low temperatures for 
an amorphous metal-nonmetal alloy system, a-Nba;Sii_a:, 
which exhibits a disorder-induced metal-insulator transi- 
tion 

Amorphous niobium-silicon alloy samples are de- 
posited on sapphire substrates by cosputtering from sep- 
arate Nb and Si sources onto rotating substrates to pro- 
duce large (19 mm diameter), thick (1 ^m) homogeneous 
samples suitable for millimeter wave transmission experi- 
ments. The sapphire substrates have their c-axis oriented 
perpendicular to the plane, and are polished so that the 
faces are parallel. At the same time other samples were 
deposited in a Hall-bar configuration for dc experiments. 
Electron microprobe analysis verifies the lateral homo- 
geneity produced by this process and was also used to 
estimate the niobium concentration (values are listed in 
Table I). 

The dc conductivity was measured for each sample in 
a standard pumped He-4 cryostat over the temperature 
range 1.4-300 K. The frequency dependent conductivity 
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is determined from millimeter-wave transmission experi- 
ments in the frequency range 100-900 GHz and temper- 
ature range 2.8-300 K. As previously described dj, the 
transmission through the sample on the substrate oscil- 
lates as a function of frequency because of standing waves 
in the substrate, which acts as a Fabry-Perot resonator 
| p5[ . We first measure the peak frequencies for each sub- 
strate alone as a function of temperature, then remea- 
sure after sample deposition. The complex conductivity 
CTi + 102 can be determined from the changes in peak 
heights and shifts in the peak frequencies p^ . 

Figure 1 shows the dc conductivity for five Nba;Sii_a; 
alloy samples plotted on a logarithmic scale versus T~^/^. 
For temperatures below 10 K, the data falls on a straight 
line, in agreement with Eq. ||. The slopes in the low tem- 
perature limit give the Tp values, which are summarized 
in Table I. Note that in the Efros-Shklovskii theory, the 
value of To becomes larger as the localization length de- 
creases; in Ref. the delocalization transition occurred 
for Nbj,Sii_:, samples with ct(77 K) > 5800 (VLm)-^ . Our 
data are not consistent with the Mott variable range hop- 
ping theory unless the effective dimension for the hopping 
electrons is cZ = 1 (which could occur if the hopping oc- 
curred preferentially along percolating paths) . 

Figures 2 and 3 show the real and imaginary parts of 
the complex conductivity respectively, for three samples 
as a function of frequency at 2.8 K. For all three samples, 
the imaginary part of the conductivity is negative (capac- 
itive) with a magnitude that is considerably larger than 
the real part. Also for all three, both real and imaginary 
parts follow a power law frequency dependence with a 
similar exponent a for both parts. This behavior implies 
that for each sample the complex conductivity can be 
described by: 
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The fact that both the real and imaginary parts individ- 
ually follow the same power law over a broad frequency 
range is compatible with the Kramers-Kronig relations 
[01 . We have three techniques for determining a from 
the data: (1) fitting 0-1(61;), (2) fitting 02 (w), and (3) fit- 
ting (72 versus cti (with frequency as an implicit variable) 
to obtain the phase angle of the complex conductivity: 



a = — tan 

TT 
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The three techniques give values for a that agree within 
experimental error, as shown in Tabic I. The value of 
a is less than unity for all three (that is, the frequency 
dependence is sublinear), and the value of a increases as 
the samples become more insulating. 

The relatively large values for the imaginary part of the 
conductivity for all three samples also implies relatively 
large values of the real part of the dielectric function 



as shown in the inset to Figure 3. For all three samples, 
the dielectric function increases slowly as frequency de- 
creases, and its zero frequency limit is not attained in our 
frequency range for Samples 1 and 2. The samples with 
higher conductivities exhibit the larger values of ei , which 
is consistent with the behavior in an insulating phase 
as a delocalization transition is approached — electrons 
in states with longer localization lengths become highly 
polarizable — as has been previously observed at radio fre- 
quencies in doped crystalline semiconductors PJl^]. 

We conclude that the electrons in these samples are 
in highly polarizable, strongly interacting, yet localized 
states. The frequency dependence of the complex con- 
ductivity possesses a nontrivial power law frequency de- 
pendence (at least in a limited frequency range that 
reaches the quantum limit). The dc conductivity fol- 
lows a temperature dependence consistent with variable 
range hopping models, so we now turn to the predic- 
tions of those models for the high frequency conductivity. 
The models can be considered to be "two site" models 
in that hopping occurs between two sites either with- 
out considering interactions at all (the Mott model) or 
without considering screening effects produced by virtual 
hops of other electrons (the E-S model). 

The two variable range hopping models produce dif- 
ferent predictions for di in the quantum regime (but not 
too high in frequency): ksT < huj < fcsTo [||J|]: 

Mott : (Ti(w) = n'^e^glhuj'^^^ \n{2Io/{riLo)) (6) 

'^~^^"^("^" 51n(2V(M) - 
Here go is the impurity band density of levels and lo is 
a microscopic energy overlap. The E-S expression as- 
sumes a broad Coulomb gap (which would be the case 
for a metal-nonmetal alloy, since the level density is quite 
high). Besides the different prediction for the frequency 
dependence (quadratic for Mott versus slightly superlin- 
ear for Efros-Shklovskii), the difference in form between 
the two formulas is striking and results from the universal 
shape of the density of levels in the Coulomb gap (which 
determines the absorbing transitions) in the Efros and 
Shklovskii theory. The Mott expression also has a much 
stronger dependence on the localization length. Note 
that the Efros-Shklovskii expression would be expected 
to depend rather weakly on the niobium concentration 
(through the logarithmic term containing the overlap in- 
tegral value) especially since the value of k is assumed 
to be determined by the host material (amorphous sili- 
con here) , which does not include the contribution of the 
hopping electrons. 

For (72, Efros |T9[ has produced a similar expression 
to Eq. I?! but without the logarithmic factor, so that (72 
should have a purely linear frequency dependence. This 
would give a frequency independent value for ei. 
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Our data for ai shows neither the quadratic frequency 
dependence nor the strong dependence of the magnitude 
of (Ji on the locahzation length predicted by the Mott 
theory is inversely proportional to Tq in the Mott the- 
ory for d = 1), which rules out attempts to explain the dc 
conductivity data as one-dimensional hopping of weakly 
interacting electrons. 

The data are not fully described by the Efros- 
Shklovskii model. We observe a sublinear frequency de- 
pendence for the complex conductivity (as opposed to 
the superlinear and linear dependences predicted for ai 
and 0-2, respectively). Also, the measured magnitudes of 
(Ti are significantly greater than the prediction of Eq. |^ 
if K is taken to be the dielectric constant of amorphous 
silicon (the dashed line in Figure 2 indicates the magni- 
tude of the numerator for k = 12; the logarithmic factor 
in the denominator is difficult to estimate but is greater 
than one in the theory). There are two caveats to con- 
sider in these comparisons, however: (1) the data may 
not be fully in the T — limit (since Hlo is only some- 
what greater than ksT for our lowest frequencies and the 
theory is for T = 0) and (2) for Sample 1, the value of 
To is low enough that the experiment does not remain 
below the upper frequency limit on the range of validity 
for Eq. over our entire frequency range. 

It may be interesting to consider the effects on the 
electrodynamics of a hopping electron in a model where 
screening effects of other hopping electrons are taken into 
account. The host material dielectric constant k used 
in the E-S model is much smaller than the contribution 
to ei from the hopping electrons themselves in our fre- 
quency range for all of our samples. Thus it is conceivable 
that the full dielectric response must be taken into ac- 
count when screening in a Coulomb glass is considered. 
Recent tunneling experiments in doped semiconductors 
have been interpreted as supporting the importance of 
many-electron composite excitations 

To conclude, we observe a power law dependence for 
the complex conductivity in the low temperature limit, 
implying that the hopping conductivity in our samples 
exhibits a richer range of behavior than can be described 
by the existing theories of variable range hopping. The 
development of theories that account for the full dielec- 
tric response of the system and additional experiments at 
lower frequencies and temperatures would appear to be 
required to fully understand the thermal and quantum 
processes that control the conductivity of the Coulomb 
glass. 
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FIG. 1. DC conductivity on a logarithmic scale versus 
T'^^'^ for a series of Nba;Sii_a; samples. Niobium concen- 
trations for the samples are given in Table I. In the low tem- 
perature limit (below 10 K) the data follow the form of Eq. 2 
(the Efros-Shklovskii variable range hopping model). 



FIG. 2. The real part of the conductivity versus frequency 
(on logarithmic scales) at a temperature of 3 K for three 
NbiSii-x samples. For all three, the data can be fit by a 

power law frequency dependence with an exponent somewhat 
less than one (the dashed line has a slope of one). Note that 
all of the measurement frequencies exceed the frequency cor- 
responding to the thermal energy: ksT/h = 60 GHz. 
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FIG. 3. The magnitude of the imaginary part of the con- 
ductivity versus frequency (on logarithmic scales) at a tem- 
perature of 3 K for three Nb^Sii-x samples. For all three, 
the sign of the imaginary part is negative (capacitive) and its 
magnitude is larger than that of the real part; the data can 
be fit by a power law frequency dependence with a power that 
agrees with that of the real part within experimental uncer- 
tainty. The inset shows the same data but expressed as the 
real part of the dielectric function, which tends to increase as 
the frequency decreases for all three samples. 



TABLE I. NbSi alloy sample parameters. The three a val- 
ues result from fitting ai and (T2 to the real and imaginary 
parts of Eq. 3, respectively, and the phase angle of the com- 
plex conductivity in Eq. 4. The value A results from fitting 
the magnitude of the complex conductivity in Eq. 4. 

Sample x a{77 K) To a&ta A 

[Nbat%] [103(fim)-i] [K] aa„a^^,a^^/^^ [(fim)-i 

±0.5 ±5% ±2% ±10% ±25% 

1 8.2 4.1 15.6 0.79,0.74,0.76 32 

2 7.8 2.9 75 0.87,0.87,0.83 16 

3 6.0 0.49 470 0.85,0.93,0.91 4.6 

4 4.9 0.176 860 

5 4.3 0.086 1010 
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